Cosmological models with fluid matter undergoing velocity diffusion 
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A new type of fluid matter model in general relativity is introduced, in which the fluid particles are 
subject to velocity diffusion without friction. In order to compensate for the energy lost by the fluid 
particles due to diffusion, a cosmological scalar field term is added to the left hand side of the Einstein 
equations. This hypothesis promotes diffusion to a new mechanism for accelerated expansion in 
cosmology. It is shown that diffusion alters not only quantitatively, but also qualitatively the global 
dynamical properties of the standard cosmological models. 
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INTRODUCTION 

The cause of several physical and biological processes 
can be attributed to diffusion. Notable examples are heat 
conduction, Brownian motion and the transport of ma- 
terials within cells. At the microscopic level, diffusion is 
due to random collisions between the particles of the sys- 
tem under study and those of a background substance. 
At the macroscopic scale, in which random effects are av- 
eraged, diffusion is described by an effective, determin- 
istic theory. The macroscopic physical quantities are so- 
lutions of partial differential equations, such as the heat 
equation or the Fokker-Planck equation. The applica- 
tions of such theories go well beyond the original phys- 
ical motivations and include fields as different as image 
processing, economics and social sciences. 

The purpose of this paper is to propose and explore 
the idea that diffusion may also play a fundamental role 
in the large scale dynamics of the matter in the universe. 
The standard assumption in cosmology is that the galax- 
ies (in which most of ordinary matter is supposed to be 
concentrated) can be represented as the particles of a 
fluid. Here it will be further assumed that the fluid par- 
ticles undergo velocity diffusion. The model will be de- 
veloped within Einstein theory of general relativity with 
a cosmological scalar field, the latter playing the role of 
the background medium in which diffusion takes place. 



THE MODEL 

The general relativistic theory of diffusion has received 
considerable attention in recent years, see [11-11]. As 
shown in [H, the energy-momentum tensor T^'^ and the 
current density of matter undergoing velocity diffu- 
sion satisfy the equations V^J^ = and 



(1) 



Here denotes the Levi Civita covariant differentiation 
on a spacetime {M,g), and a is a positive dimensional 
constant. In ref. [H the matter field is given by the one- 
particle kinetic distribution function / in phase space, 



and the tensor fields T^'^, are defined by suitable in- 
tegrals of / in the momentum variable. The relation ([T]) 
is a consequence of the evolution equation satisfied by 
/, which is a degenerate Fokker-Planck equation without 
friction on the tangent bundle of spacetime. Hence the 
model studied in |1] describes the kinetic motion of an 
ensemble of particles undergoing diffusion in the velocity 
variable. The right hand side of ([TJ is a macroscopic rep- 
resentation of the diffusion forces acting on the particles, 
and (7 is the diffusion constant. 

This paper is concerned with fluid matter models with 
energy-momentum tensor and current density 



that satisfy 



(2) 



(3) 
(4) 



The scalar functions p,p,n represent respectively the 
rest-frame energy density, the pressure and the number 
density of the fluid, while the vector field : u^u^^ = — 1, 
is the four-velocity of the fluid. Motivated by the kinetic 
model introduced in [ij, an ideal fluid that satisfies ([1]) 
will be said to undergo (molecular) velocity diffusion and 
the parameter a will be called diffusion constant in the 
present context as well. 

Eq. ^ is equivalent to the following equations 



V^(pu^ 



an. 



{p + p)u^'\7^u'' + u'^u^'Vf.p + V> = 0, 



(5) 
(6) 



which are obtained by projecting ([3]) into the direction of 
and onto the hypersurface orthogonal to u'^, respec- 
tively. Note that the Euler equation ^ is the same as in 
the diffusion-free case; only the continuity equation ([5]) 
is affected by the presence of diffusion. This is due to the 
fact that the diffusion force anu^ acts on the direction 
of the matter flow. 

To transform ([4])-(l6l) into a complete system on the 
variables (p, p, n, u^), an equation of state must be added. 
Letting the energy per particle e and the specific volume 



2 



V be defined as e = p/n, and v = 1/n, the first law of 
tfiermodynamics states that 



de — —pdv + Tds, 



(7) 



where the scalar functions T and s denote respectively 
the temperature and the entropy of the fluid. Units are 
chosen so that Boltzmann's constant is equal to one. An 
equation of state relating the thermodynamic variables 
may be prescribed in the form 



p = h{n,s), 



(8) 



for some function h, which for simplicity will be assumed 
to be continuously diffcrcntiable. Substituting in ([7]) we 
obtain 



p + p 



dh 
dn ' 



T(n, s) 



1 dh 
n ds 



(9) 



For ordinary fluids the specific enthalpy c and the tem- 
perature T are required to be positive, which implies that 
p + p > 0, dh/dn > 0, and dh/ds > 0. Then by (0) 
and (jlj we obtain 



u^V.s — -. 

^ T(n,s) 



(10) 



The equations ©, ©, dH), ®, ^ provide the desired 
complete system for describing the evolution of the mat- 
ter field variables {n,s,u^). Eq. ([5]) can be replaced by 
eq. ([ll, since they are equivalent when (fTO|) holds. Note 
also that the entropy of the fluid is strictly increasing 
along the matter flow, unless diffusion is switched off (i.e., 
a is set to zero). This is due to the fact that diffusion is 
an irreversible process. 

For the applications in cosmology, the case when the 
pressure p and the rest-frame energy density p are con- 
nected by a liner relation p = (7 — l)p, 7 = const., is 
particularly important. By ^ this can be achieved by 
choosing 



h{n, s) = nTls), 



(11) 



where 7 > to ensure e > and F'{s) > to ensure 
positivity of T. In this situation, the entropy variable 
s can be replaced by the modified entropy S = F{s) 
and (jlOp simplifies to 



(12) 



The consistent description of any matter dynamics in 
general relativity requires the coupling with the Einstein 
equation: 

R^.. ' ^g^^.R = Tf,, {8ttG = c=1), (13) 

where i?^^ is the Ricci tensor of the metric g and R = 
g^'^Rfj.u- However it is clear that fiuid matter undergo- 
ing velocity diffusion cannot appear as the only source 



in the Einstein equation. In fact the Bianchi identity 
^^{Rfii/ — ^gtj.uR) = implies that solutions of the Ein- 
stein equation can only exist if V^T^"^ — 0. The phys- 
ical reason for this incompatibility is that our new fluid 
model does not take into account the interaction between 
the particles undergoing diffusion and the background 
medium in which diffusion takes place. While the dif- 
fusion approximation in Newtonian mechanics consists 
exactly in neglecting this interaction [H, we see that in 
general relativity such an approximation, which leads in- 
evitably to a model that does not preserve energy, is in- 
consistent. The next best "diffusion approximation" con- 
sists in assuming the simplest possible interaction with 
the background medium. An interesting choice, because 
of its simplicity and because it leads to a natural gener- 
alization of the cosmological constant assumption, is to 
model the effect of the background medium by a scalar 
field according to the following modification of (fT3)) : 



Rfj.i' 2 gt^fR 



= T^, (8^G = c=l). (14) 



The scalar field (j) will be called cosmological scalar field, 
since it enters in (ITil) in the place where the cosmolog- 
ical constant A is commonly introduced. It should be 
emphasized once again that the cosmological scalar field 
is not to be interpreted as ordinary matter, but rather 
as a background field which, interacting with the fluid 
particles, causes their diffusion. According to particle 
physics, 4> plays the role of the vacuum energy. 

By Q , the cosmological scalar field satisfies the equa- 
tion 



anu^ 



(15) 



There are a number of important properties that one can 
infer from (fTS]) . Firstly, (fT5|) forces the matter flow to be 
irrotational: = V ^^{nu^) — \/^(nu^) = 0, where uj^i, 
is the vorticity tensor of the fiuid. Secondly, <j) satis- 
fies the transport equation u^SJ ^cj) = —an, which shows 
that (f) is decreasing along the matter flow. Finally, since 
the fluid number density is assumed to be a conserved 
quantity, i.e., ^{nu^) = 0, the cosmological scalar field 
satisfies the homogeneous wave equation V^V^(/) = 0. It 
follows that the canonical energy-momentum tensor of 
the scalar field is divergence- free: 



Note however that this energy-momentum tensor is not 
added as a source in the Einstein equation (ITU) . This 
again is due to our interpretation of as a background 
medium, rather than an ordinary matter field. The en- 
ergy lost by the particles due to diffusion is stored in 
the scalar field </> and is propagated through spacetime 
in form of waves without dissipation. In the absence of 
diffusion, i.e., when a = Q, (pis constant and (1141) reduces 
to the Einstein equation with cosmological constant. 
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COSMOLOGICAL SOLUTIONS 

Next we consider spatially homogeneous and isotropic 
solutions, for which the line element can be written in 
the form 



ds'^ = -dt^ + a{tf 



dr^ 



1 — kr^ 



r'^dn^ 



(16) 



where a{t) > is the scale factor, fc = or ±1 is the 
curvature parameter of the hypersurfaces t = const, and 
drt is the invariant surface element on the unit sphere. 
This class of solutions is the most important and widely 
used for the applications in cosmology @. Let us write 
the fluid equations derived previously on a Lorentzian 
manifold with line element ([T6l) . In this case the matter 
field variables are functions of t only. Moreover we as- 
sume that p = h{n, s) is given by (ITT|) . which implies that 
p = (7 — for some constant 7. Although this is not 
always necessary in what follows, it will be assumed that 
2/3<7<2;in particular the fluid satisfies the strong 
and the dominant energy condition. The conservation of 
particles number, eq. (|3|), implies 



(17) 



where a subscript means evaluation at time t = 0. The 
equation (|12|) on the modified entropy S becomes 



1-7 



a{t) 



3-37 



(18) 



where a dot denotes differentiation with respect to t. 
Next the equations on the geometric fields will be consid- 
ered. The equation (jisp on the cosmological scalar field 
6 becomes 



while the Einstein equation (|14p entails 



H 



k 



a{tf 



>7 - 1 P 



where 



a 
a 



(19) 

(20) 
(21) 

(22) 



H denotes the Hubble function. For any specific choice of 
7 and fc = 0, ±1, the equations (fT7|) - ([2T|) form a complete 
system for the fluid variables {n,S), the scale factor a, 
the Hubble function H and the cosmological scalar field 
(j). The initial data consist of ag > 0, no, So,(t>o, Hq such 



that pOj) is satisfied at time t = 0. Identifying t = with 
the present value of the cosmological time, it is physically 
justified to assume Hq > — since the universe is cur- 
rently expanding — and (po > — since (f)o corresponds 
to the present observed value of the cosmological con- 
stant. It will also be assumed of course that no, So > 0. 

In the following the global behavior of solutions 
to pT|) - (PT|) in the flat case (fc = 0) will be analyzed. Let 
us start by recalling the global behavior of the analogous 
diffusion- free cosmological model (cr = 0). In the ab- 
sence of diffusion, one recovers the Einstein-perfect fluid 
system with non- negative cosmological constant A ~ (j)Q. 
In this case it is well known that the model develops a 
singularity in flnite time in the past, while in the future 
it is singularity-free and forever expanding. (In fact, one 
can even solve explicitly the equations in this case, see [^, 
Sec. 14.2].) It will now be shown that when diffusion is 
switched on, the past and future qualitative behavior of 
the universe may be completely different. 

Since fc = 0, we can introduce cartesian coordinates 
{x,y,z) such that ds^ = —dt^ + a{t)'^{dx^ + dy"^ + dz^) 
and oq = 1. 

Since 0o > and (f) is decreasing, the cosmological 
scalar held is always positive in the past, causing ac- 
celerated expansion of the universe, and could become 
negative at later times. If this happens the universe is 
doomed to collapse into a singularity in finite time. In 
fact if (f>{t) < at some time t, then by (PT|) we will have 
H < —\(f){t)\/3'-H'^ for all t > f and therefore there exists 
a time t* such that _ff (i) — )■ — 00 as t ^ i*, which implies 
a(t) ^ as f — > i*. On the other hand, if (f>{t) never 
vanishes, then the cosmological model will be singularity 
free and forever expanding in the future. This is a trivial 
consequence of l|^D|) t.-n. which shows that H never van- 
ishes if > 0, and therefore, since Hq > 0, it will remain 
positive for all times. Let (jjaa = limt_>oo > 0. It will 
now be shown that when (j)oo > 0, the metric behaves 
asymptotically as the De Sitter solution, precisely 



a[t) exp \ — — 




as t 



00. 



(23) 



In fact, by L'Hopital's rule. 



S 



S 



a 



a 



Moreover by ((20)) - (|2T|) . 
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Whence a^H — ^ 00, and thus p ^ 0, as t ^ 00. By 
by (Uni), ^ 0oo/3, which implies (|^ . 

The behavior of the cosmological model in the past is 
more sensitive to the choice of the initial value Sq rather 
than 00. If So is sufficiently small, the modified entropy 
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(a)Futurc collapse (Sq = 3,4>o = 0.3) 
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(a)Futurc collapse (5o = 0.6, (/)o = 0.1) 
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(b)Unlimited expansion (So = 3, 0o = 1) 
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(b)Unlimited expansion (So = 0-6, 4>q = 1) 



FIG. 1: Examples of solutions with big-bang singularity (a 
7 = no = ao = 1) 



FIG. 2; Examples of solutions without big-bang singularity 
(a = 7 = no = ao = 1). The solutions can be matched to De 



Sitter at the time t — to 



-0.6. 



S could vanish at some negative time tg, while a{to) is 
still positive. (For instance, in the case of dust (7 = 1), 
the solution of ([T^ is S{t) = So + crt and therefore to = 
—So/c.) In this case the solution is unphysical for t < to, 
because p would be negative and the entropy decreasing, 
even if the metric of spacetime remains smooth. However 
a physically viable and interesting cosmological model 
can be constructed by matching the metric at the time to 
with the De Sitter solution aosit) = Ccxp(-\/ (j){to)/3t), 
where the constant C is such that aDs{to) = a{to). The 
resulting cosmological model has no big-bang singularity 
and is vacuum up to the time to (since p = and a{t) = 
o.Dsit) for t < to), at which time the vacuum energy 
^(io) starts to be converted into matter energy p. The 
other possibility is that the scalar factor vanishes while 
the modified entropy is still positive. In this case a big- 
bang singularity forms in the past. 

Figures [1] and [2] show examples of solutions obtained 
numerically. The solutions depicted in Figure 1 have a 
big-bang singularity in the past, while those on Figure 2 
do not. Moreover, the solutions depicted are all dust so- 
lutions (7 = 1) and the constants a, no have been fixed to 
one. Thus S{t) = So + 1 in all pictures. Similar solutions 
can be obtained for fc = ±1 and for all other values of 
7 G [2/3,2) and cr > (no matter how small). Whether 
any of these solutions may represent the actual behavior 



of our universe requires of course experimental evidence. 
In particular, while observational constraints are avail- 
able for the values of (no, po,4>o. Ho), the possibility that 
cr > remains open. 

We conclude by pointing out that a simple explicit so- 
lution of (fT7 l) - ([2T1) for fc = and 7 = 1 is given by 



a{t) = 1 + (-|-) t, 
(nool2fl^ 



(24a) 



no 



2\ 1/3 



a{t). (24b) 



The only free parameter in the solution (p4)) is no, 
thus it is highly non-generic. The qualitative behavior 
of this solution is similar to the one depicted in Fig- 
ure IHb) . Note also that for cr = the solution reduces 
to {a,(j),S) = (1,0,0), which corresponds to Minkowski 
spacetime. Thus the solution (IM)) is a "pure diffusion" 
one. 
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